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Abstract 

We test several BFKL-like evolution equations for unintegrated gluon distributions against 
forward-central dijet production at LHC. Our study is based on fitting the evolution scenarios 
to the LHC data using the high energy factorization approach. Thus, as a by-product, we 
obtain a set of LHC-motivated unintegrated gluon distributions ready to use. We utilize this 
application by calculating azimuthal decorrelations for forward-central dijet production and 
compare with existing data. 


1 Introduction 

A typical procedure in applying QCD to hadronic collisions relies on factorization theorems. They 
consist in two ingredients: a perturbatively calculable hard part and a nonperturbative piece 
parametrizing hadrons participating in a collision. The most known and tested is the collinear 
factorization (see e.g. ^ for a review), which applies for a variety of processes, including jet 
observables in deep inelastic scattering (DIS) and hadron-hadron collisions. Here, the nonper¬ 
turbative component is parametrized in terms of parton distribution functions (PDFs) which un¬ 
dergo Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations. The key feature 
of PDFs is the universality, i.e. the PDFs that are measured in one process can be used in any other 
for which the factorization holds. Therefore, for instance one can use PDFs fitted to DIS structure 
functions and use them to make predictions for jets in hadron-hadron collisions. Although the 
collinear factorization is powerful and well-tested, it is supposed that for certain observables, e.g. 
forward jets at high energies, another kind of evolution equations for the PDFs is needed. Namely, 
the perturbative calculations contain the logarithms of the form ag log(l/a:), where x is the lon¬ 
gitudinal fraction of the hadron momentum carried by the parton. At high energies and forward 
rapidities x is small and these logarithms need to be resumed. This is accomplished by means of 
various “small x” evolution equations, which essentially are various extensions of the pioneering 
Balitski-Fadin-Kuraev-Lipatov (BFKL) evolution equation (see e.g. [2]). In the small x domain 
the transverse momenta of the partons exchanged between the perturbative and nonperturbative 
parts are not suppressed comparing to the collinear factorization. Therefore, the PDFs have an 
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explicit dependence on the transverse momentum of a parton. Such objects are often referred to, 
as transverse momentum dependent PDFs (TMDs) or Unintegrated PDFs, although the former are 
typically used outside the small x physics, and posses unambiguous (though in general process de¬ 
pendent) field theoretic definitions. Actually, at small x one usually deals with initial state gluons 
only, and thus the object of interest in this paper is an Unintegrated Gluon Distribution (UGD). 
The UGDs have to be convoluted with a perturbative “hard part” according to so-called kx or High 
Energy Factorization (HEF). We describe this approach in some more detail in Section]^ Here, let 
us just mention that unlike the collinear factorization, the HEF is not a QGD theorem and actually 
the universality of UGDs is supposed to be violated for jet production in hadron-hadron collisions. 
Thus, in principle, the standard procedure of fitting the UGDs to the F 2 HERA data and using 
it for jets in hadron-hadron collisions is not correct, but there are no quantitative measures of the 
factorization violation so far. Actually, HEF is surprisingly quite successful with describing LHG 
data using UGDs from fits to structure functions, see for instance [3]. At present, there are several 
fits to F 2 data using different small x approaches, see |1H7] for more details. 

In the present work we undertake another path. We make an attempt to fit various BFKL- 
like UGDs directly to the LHG data for jet forward jet production. It has a twofold purpose. 
First, we have an opportunity to explore UGDs using relatively exclusive observables. Second, 
we want to free ourselves from the aforementioned universality problem when transferring UGDs 
from DIS to the LHG domain. We consider two separate measurements: jet transverse momentum 
spectra [S] in forward-central jet production and forward-central dijet decorrelations j^. The first 
measurement consists of two separate sets of data: for the forward jet and for the central jet. Thus, 
the mutual description of both spectra imposes a strong constraint on the UGDs and we shall use 
this measurements to make our fits. The second measurement will be used to test the fits. 

The paper is organized as follows. In Section]^ we describe the approach of HEF. The small 
X evolution equations with various components incorporating sub-leading effects are discussed in 
Section The fitting procedure and the software used are described in Section We give the 
results in Sectionj^ Having the fits, we test them against recent forward-central dijet decorrelations 
data in SectionFinally, we discuss our research in Section]^ 


2 High Energy Factorization 

In this introductory section we discuss in more detail issues concerning factorization at small x. 
This task is somewhat complicated, notably because of the various existing approaches and various 
existing definitions of UGDs. 

In the following paper the notion of HEF corresponds to a general class of factorization ap¬ 
proaches supposed to be valid at small x. Below we list some of the existing realizations: 

1. the factorization of Gribov, Levin and Ryskin (GLR) [TU] for high-pr inclusive gluon pro¬ 
duction 

2. the factorization of Gatani, Giafaloni and Hautmann (GGH) |llllT3] for heavy quark produc¬ 
tion in DIS, photo-production and hadron-hadron collisions 

3. the factorization of Gollins and Ellis m for heavy quark production in hadron-hadron col¬ 
lisions 

4. the factorization for inclusive gluon production in the saturation regime for proton-nuclei 
collisions within the Golor Glass Gondensate (GGG) approach [14] and color dipole formalism 
|T51 [T3] (the equivalence of both approaches was shown in ini) 
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In these approaches the nonperturbative part is parametrized in terms of UGDs undergoing BFKL 
evolution (for GRL, GGH, Gollins-Ellis factorizations) or nonlinear Balitsky-Kovchegov evolution 
[m [19] (for GGG). On the other hand, superficially similar objects to UGDs appear in so-called 
transverse momentum dependent (TMD) factorization and are called TMD PDFs. One should 
however realize that the enumerated approaches are valid at leading logarithmic approximation, 
while the TMD factorizations are valid to all orders in the leading twist approximation. Moreover, 
unlike most of UGDs in the HEF factorizations, the TMD PDFs have precise operator definitions in 
terms of matrix elements of nonlocal operators. Those definitions require appropriate Wilson lines 
to be inserted in order to make the definitions gauge invariant and to resum collinear gluons related 
to final and initial state interactions. These insertions make the TMD PDFs, in general, process 
dependent and thus non-universal, breaking the principle of factorization (for more details see e.g. 
EOlllI])- Only for processes with at most two hadrons the TMD factorization is proved to hold to 
all orders (for example back-to-back single hadron production in DIS or Drell-Yan scattering). The 
natural question arises whether the non-universality of TMD PDFs transfers to the small x limit. 
In ref. (Hj an explicit arguments were given that this is the case for dilute-dense collisions (actually 
the arguments hold for so-called “hybrid” factorization - see also below). Moreover it is known 
from the GGG approach that at really small x, i.e. in the saturation regime, the cross sections 
cannot be described by just dipoles (averages of two Wilson lines), but also higher correlators 
are needed [23], what violates the ordinary logic of factorization. However, for the case of back- 
to-back dijet production in dilute-dense collisions a generalized factorization has been proposed 
|24| ; that is, the cross section can be given in terms of hard factors and certain universal pieces. 
Recently, these results were improved to the case of imbalanced dijets |25j . In particular, when 
the imbalanced transverse momentum is of the order of transverse momenta of the jets the HEF 
for dijet production can be derived from the dilute limit of the GGG approach. 

In the present work we shall constrain ourselves to dijet production in p-p collisions in the 
linear regime, as the kinematics we are interested in (and where the data exist) do not allow to 
develop the saturation region. We want to utilize most of the phase space covered by the data, 
thus we do not constrain ourselves to the back-to-back dijet region analyzed in m- Rather, we 
shall use the HEF factorization for dijet production. Since this approach is an extension of the 
GGH formalism, we shall now briefly recall the latter and the required extensions to obtain HEF 
for dijets. For a direct derivation from GGG approach see |25| . 

In the GGH high energy factorization, one considers the heavy quark pair produced via the 
tree-level hard sub-process g* {kA) g* {ks) QQ in the axial gauge. The initial state gluons are 
off-shell and have the momenta of the form = xaPa + kxA and ks = XbPb + kxB, where pA, 
Pb are the momenta of the incoming hadrons and pA ■ kxA = Pb • kxB = 0. This particular form 
of the exchanged momenta is a result of the imposed high energy limit. The off-shell gluons have 
“polarization vectors” that are pA and pB respectively. Thanks to this kinematics the sub-process 
given by ordinary Feynman diagrams is gauge invariant despite its off-shellness. In GGH approach 
the factorization formula for heavy quark production reads (see Fig. S') 

i^AB^QQ = II ^ I XkrB j ^ 

J^g*/A{xA,kTA) J^g*/B ixB,kTB) da^,( xa, Xb, kxA, kxB) , (1) 

where dag,g ,^qq is the partonic cross section build up from the gauge invariant g*g* -A QQ 
amplitude and Tg*/ a, d^g-/B are UGDs for hadrons A and B. The contributions with off-shell 
quarks are suppressed. The UGDs are assumed to undergo the BFKL evolution equations. In 
Ref. [12] it was argued that similar factorization holds to all orders for DIS heavy quark struc¬ 
ture function, although the argumentation misses the details comparing to collinear factorization 
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Figure 1: A) The CCH factorization for inclusive heavy quark production; despite the fact that 
the gluons entering the central blob are off-shell the sub-process is gauge invariant. B) For sub¬ 
processes with final state gluons the gauge invariance requires the off-shell gluons to be replaced 
by the effective particles giving rise to multiple eikonal gluon exchanges between the blobs. 


proofs [T], especially the definitions of UGDs and complications arising at higher orders in the axial 
gauge pS] . 

In the works as well as in this paper the CCH factorization was extended 

to model the cross section for jet production in hadron-hadron collisions. The first difficulty arises 
because now one has to consider also gluons in the final state, e.g. g*g* —)■ gg sub-process for 
dijet production. The corresponding amplitude is however not gauge invariant when calculated 
from ordinary Feynman diagrams. A few approaches have been proposed to calculate a gauge 
invariant extension of such amplitudes These gauge invariant off-shell amplitudes in 

fact correspond to a vertex that can be calculated from the well-known Lipatov’s effective action 
|35l [M| (see Fig. [^). The approaches IMIf34] were however oriented on practical and efficient 
computations of multi-particle off-shell amplitudes using helicity method and computer codes. As 
stated before, in CCH the UGDs undergo BFKL evolution. In our extensions of CCH approach 
we allow the UGDs to undergo more complicated evolution equations, which are more suitable for 
jets. More details will be given in Section Yet another modification of the CCH formula comes 
from the fact that the present study concerns the system of dijets where one of the jet is forward, 
while the second is in the central region. From 2 —>■ 2 kinematics it follows then, that xa ^ xb (or 
the opposite), except for the small corner of the phase space. Since is typically of the order of 
0.5 the usage of small x evolution for iFg */b is questionable (this is similar to dilute-dense system 
considered e.g. in [H]). Therefore we use collinear approach on the B hadron side Technically, 
one takes the collinear limit in dag-g*^ 2 j by sending kxB —^ 0 to obtain a sub-process with one 
off-shell gluon dag*g^ 2 j (the off-shell amplitudes have well defined on-shell limit). In this limit 
one has to take into account also sub-processes with initial state on-shell quarks, dcFg*q^ 2 j- The 
remaining integral over d^kxB gives helicity sum for B partons on one hand, and the integrated 
(collinear) PDF on the other J dk'^ {xBjkxB) = fa{xB)- Thus, the final formula for the 

factorization model reads 



^q* [A i^XA , kxA ^ /^) fb {xb ^ /^) dCg^fj — y2j (^XA , Xb , kxA ; fd) ^ (2) 

where we have included the hard scale dependence not only in the collinear PDFs fb, but in the 
UGD as well. Such a dependence turns out to be important for certain exclusive observables 
involving a hard scale (e.g. large px of jets; see e.g. El)- We note, that when the final states 
become well separated in rapidity, i.e. when the central jet lies in the opposite hemisphere to the 
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forward jet we start to violate our condition xa Xb and different approach should be used. 
The factorization formula ([^ resembles the linearized approach of [23] but it extends beyond the 
correlation limit as here the hard sub-processes have injected a nonzero kx- As mentioned before, 
the formula ([^ has been recently derived from the CGC approach in |25| . 


3 Small X evolution equations 


Let us now discuss the evolution equations for UGDs which were used in our fits. As described in 
the preceding section we concentrate on linear evolution equations. Below we list some of them 
with a short explanation. We consider only gluon UGDs, thus we skip the subscripts in Xg* /a- 

1. pure BFKL equation 

The equation in the leading logarithmic approximation reads |38l I39| 


F (a;, kx) = Fq (a;, kx) 

(f ’ 9 t) ~ (f ’ ^r) , (f > ^t) /on 
\qT-kx\ V^qT + K . 

where Og = Ncag/'x with W being the number of colors. The initial condition for the 
evolution is given by Fq. The NLO BFKL equation is also known gsi El. One of the 
drawbacks of the pure BFKL equation comes from the fact that qj. of the gluons emitted along 
the ladder is unconstrained. Indeed, since in the BFKL regime the virtuality of the exchanged 
gluons is dominated by the transverse components, the resulting kinematic constraint reads 

131113] 

Qx < -—- kx ^ - kx- (4) 

z z 

This constraint is also often referred to as the consistency constraint. 

2. BFKL with the kinematic constraint (BFKL+G) 

To incorporate the consistency constraint one may include the appropriate step function into 
the real emission part of the BFKL. This operation, actually introduces some higher order 
corrections into the BFKL equation [33] ■ In addition, one may introduce another class of 
sub-leading corrections by allowing the strong coupling constant to run with the local scale 
along the ladder. Finally, one may define the q'^ integration region to lie away from the 

infrared nonperturbative region by separating the dq^ integration and moving it to the 
initial condition (the infrared cutoff k^g is taken to be of the order of 1 GeV). The improved 
equation reads |44] 



F(x,kx) — Fq(x, kx) 

qx^ (f Q (^T ~ ^gr) ~ (f ^ 

\qT-kx\ V^qT + kx 

(5) 

Recently, it has been studied in the context of the Mueller-Navelet jets, that the energy- 
momentum conservation violation (which above is cured by a “brute force”) becomes less 
harmful when full NLO corrections are applied [35] . The effects of the kinematic constraints 
in the approximate form Q as well as in the full form have been recently analyzed |46| in 
the context of the GGFM evolution equation [471 - 150] . 
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3. BFKL with the kinematic constraint in re-summed form (BFKL+CR) 
The equation © can be casted in yet another form m 


T (x, fcy) = ^0 (x, k^) 


+ (Xg (^ t ) / — 


dz f°° d^qr , 2 2\ a ^ 1.2 2 \ t I ^ 

Z I u2 TTQn-i 
'X J 


kr + qt 


where 


Ar (z, k^, fi^) = exp In ^ In 


( 6 ) 

(7) 


is the so-called Regge form factor. This form has been used in Ref. m to propose a 
non-linear extension of the CCFM equation. The scale /r has been introduced to separate 
unresolved and resolved emissions in ([^, i.e. the emissions with and q^ > fp, and 

further the unresolved part was re-summed to obtain the Regge form factor. Note, that the 
UGDs undergoing this equation do not explicitly depend on the scale ^ and that the new 
form of the initial condition has to be used (this is denoted by a tilde sign). 


4. BFKL with the kinematic constraint and DGLAP correction (BFKL+GD) 

In Ref. IH] yet another improvement of ^ was proposed. One can make an attempt to 
account for DGLAP-like behaviour by including the non-singular part of the gluon splitting 
function (the third term below) 


(x, k^) = Fq (x, k^) 

(f ’ dr) ® (^r ~ ^9 t) ~ k^F (|, q^) k^F (|, k^) 
kr-ferl \/Mt + K _ 

+ as ( 4 ) ^ ^ (^^99 (^) - 1 ) dqlF (^, 4 ) , ( 8 ) 



where Pgg (z) is the standard gluon splitting function. This correction, similar to the kine¬ 
matic constraint, accounts for certain sub-leading corrections to the BFKL equation. 


5. BFKL with DGLAP correction alone 

This variant is used to test the significance of the DGLAP term alone. 


The above UGDs do not involve any hard scale dependence. For observables involving high-pr jets 
a presence of large scale ^ in perturbative calculations would involve additional logarithms 
of the type log which can spoil the procedure. Therefore a re-summation of those logs is 

desired and it accounts in hard scale dependence for UGDs, c.f. Eq. ([^. The approach which 
incorporates both x, k^ and p? dependence in UGDs is provided for example by the GGFM 
evolution equation (the code available for a practical use is described for example in |S2])- Another 
approach, so called KMR (Kimber-Martin-Ryskin) procedure |53l [54] . takes ordinary PDFs and 
injects kx dependence via the Sudakov form factor taking care of matching to the BFKL evolution 
at small x. A serious advantage of this procedure is that one can use well known PDF sets, fitted to 
large data sets. Yet another approach was used in in therms of so-called “Sudakov resummation 
model”. This procedure reverts, in a sense, the logic used in the KMR and uses the Sudakov form 
factor to inject the hard scale dependence instead of kx- The procedure is parton-shower-like, i.e. 
it is applied after the MG events are generated and the cross section is known, and is unitary (i.e. 
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the procedure does not change the total cross section). The advantage is that one may use it on 
the top of UGDs involving nonlinear effects. The basic idea behind the model is that it assigns the 
Sudakov probability P for events with given kx and a hard scale ^ ^ pr- Then, the probability of 
surviving is 1 — P. For events with small kx and large p the emission probability is P ^ 1 and the 
unitarity of the procedure transfers such events to the region kx px- There is one more approach 
proposed in Ref. |5S], similar to the one just described, where analogous procedure is applied at 
the level of UGDs by fixing its integral over kx (it has an advantage of being independent on any 
software and one may produce grids for a practical usage). In summary, we may consider the 
following modifications of UGDs m 

6. BFKL with the Sudakov (BFKL+S) 

7. BFKL with the kinematic constraint and the Sudakov (BFKL+GS) 

8. BFKL with DGLAP correction and the Sudakov (BFKL+DS) 

9. BFKL with the kinematic constraint in re-summed form and the Sudakov (BFKL+GRS) 

10. BFKL with the kinematic constraint, DGLAP correction and the Sudakov (BFKL+GDS) 

Unfortunately, as far as fitting of UGDs is considered, the above Sudakov-based models are not 
suitable. This is because they require the knowledge of an integral (whether it is a cross section 
or integrated gluon, c.f. |3] vs [55]) which is unknown at the stage of fitting. In principle, one 
could try to use the method of successive approximations with the Sudakov model of Ref. |3]- We 
shall report on our attempts in Section There is one more comment in order here. The Sudakov 
resummation model is very sensitive to the region < 1 GeV which is not well described by the 
practical implementations of the equations [l][^ as they use certain Iow-^t cut, kxo- For kx < kxo 
the UGD is typically modelled or extrapolated by a constant value. 

Let us now discuss the models for the initial condition Pq. In this paper we have tested the 
following models (in the brackets we give the aliases used below to identify the model): 

A. exponential model (EXP) 

Po {x, kx) = N (1 — a;)“ (1 — Dx) (9a) 

B. (negative) power-like model with running Us (POW) 

Po {x, kl) = Nx^{l- xY (1 - Dx) (9b) 

fXrp 

C. DGLAP-based model (Pgg) 

Po (x, k^) = J dz Pgg {z) Go {x) , (9c) 

where 

^0 {x) = N x^ {1 — x)°' (1 — Dx) (9d) 

is a model for an integrated gluon density. 

The parameters N, A, a, D are, in general, free parameters and need to be fitted. 

We see that in principle there are quite a few variants to be fitted. Though not all of the 
combinations make sense, we are still left with several scenarios to be tested. 


7 




4 Fitting procedure 


We have used two data samples measured by CMS detector jS] for inclusive forward-central dijet 
production at CM energy ^/s = 7TeV. The central jet is defined to lie within the pseudo-rapidity 
interval |? 7 c| < 2.8 while the forward has to lie within 4.9 > \r]f \ > 3.2. Both jets are high-px jets 
with px > 35 GeV. The jets were reconstructed using anti-fcT’ algorithm with radius R = 0.5. The 
data samples consist in jet px spectra for forward and for central jets, das / dpx^Tjs with S = f,c. 
There are in total 12 data bins for both forward and central jets. 

We have applied the following fitting procedure. For each existing experimental data bin B 
we produce a 2-dimensional normalized histogram with bins in x and kx, such that the cross 
section can be calculated as 

where i,j enumerate the bins in (x, kx). To make the histograms we 


1. generate Monte Carlo events for the process under consideration with F = R*, where F* is 
a relatively “broad” trial UGD (evolving according to one of the scenarios Eng. 


2. make histograms in [x, kx) of contributions to each data bin B, 

3. divide by F* {x, kx), i.e. = i)fj/F* {x (i ), kx (j)). 

Hence, in principle R^ are independent of F* used for their generation and are calculated only 
once. This is advantageous, as the hard cross section calculation is costly in CPU time. The latter 
is calculated using the Monte Carlo C++ program Lx Jet |H5] implementing The generated 
events (weighted or unweighted) are stored in a ROOT file for further processing. For the 
UGD evolution according to scenarios [T][^ we solve the corresponding integral equations by a 
straightforward numerical iteration over a grid over x and kx- 


In order to make the fitting feasible, we need a fast routine to calculate F used in (101 for 
the cross section calculation. However, since our numerical procedure is too slow for that, we 
prepare grids over which we can interpolate the fitting parameters. Each such grid corresponds to 
a particular parametrization model and arguments range. Out of four parameters (iV, A, a, D) 
of the initial conditions, we fix D = 0 (see Sec. [^. Moreover, we note that the solution for F is 
linear in N. Thus the actual grids are in A and a. 


5 Results 


We have applied the procedure described in the preceding section to most of the models m and 
initial conditions A]|C The best values of x^/NDP (y^ per data point) are listed in Table |l| for 

mi 


models 


Note, that some of the scenarios were unable to describe the data, in particular the 


pure BFKL and BFKL with the kinematic constraint only. Evidently, the DGLAP correction is 
essential. The fitted values of the parameters of the initial conditions, N, A, a, for scenarios with 
X^/NDP < 2 are collected in Table The fits are presented in Figs. |2][^ For a better comparison 
we also plot the cross-sections scaled by p^. We observe that all the models with the DGLAP 
correction give excellent description of the central-jet data, while the px spectrum of forward jets 
is reasonably reproduced though less accurately. We also note that the models with lowest 
result in very similar predictions for the px spectra. 

Our attempts to fit the scenarios with the Sudakov resummation can be summarized as follows. 
First, we observe that the model has a small overall effect on the px spectra, although it slightly 
shifts the theory points away from the data points. We illustrate this in Fig. where we applied 








^0 

BFKL 

BFKL+G 

BFKL+D 

BFKL+GD 

BFKL+GR 

EXP 

2.4 

2.2 

1.24 

1.11 

1.52 

POW 

2.3 

1.9 

1.02 

1.12 


Pgg 

- 

- 

1.13 

1.11 



Table 1: The values of x^/NDP for fits of unintegrated gluon density evolving according to various 
models described in Section]^ The first column lists the initial condition ansatz, see also Section 
for details. 


model 

N 

A 

a 

BFKL+GR (EXP) 

0.095 

0.012 

0* 

BFKL+D (EXP) 

0.37 

0.18 

0.5* 

BFKL+GD (EXP) 

0.68 

0.14 

2.5* 

BFKL+G (POW) 

320 

1.4 

61.0 

BFKL+D (POW) 

12.7 

0.5* 

5.7 

BFKL+GD (POW) 

562 

0.96 

35.7 

BFKL+D (Pgg) 

106 

1.2 

2.5 

BFKL+GD (Pgg) 

628 

2.9 

5.7 


Table 2: The values of initial condition [Al|C] parameters obtained from the fits to the CMS data. 
We list only the scenarios with y^/NDP < 2. The values denoted by a star were fixed — see the 
main text for details. 


the Sudakov model on the top of the events obtained with one of the fits. When we now try to 
refit the parameters, we change the total cross section (used already to apply the resummation) 
and the fit fails. Although we observe that the successive iterations improve the fit, the procedure 
turns out to be insufficient to make a reliable fit with the Sudakov resummation. 

A few comments are in order. The considered jet data are not sufficient to precisely determine 
all the parameters {N, A, a, D) of the initial parametrizations ([^. Thus, first we neglect the 
(1 — Dx) factor, i.e. we take D = {). We have checked that we get no improvement when D 
is a free parameter. Next, in some cases the fits are not sensitive enough to uniquely determine 
the three remaining free parameters. In these cases we fix A or a at some plausible value (these 
are marked with a star in Table |^. Actually, besides the initial condition parameters N, A, 
a, D we have also the boundary values of kinematic parameters xa, (c.f. Q), which - to 
certain extent - are free parameters as well. We set them as follows. First, in order to be in 
an accordance with the assumptions leading to (§ we imply the cut xa < Xb- Next, for all 
scenarios we set XAmin = 0.0001. For the model with the DGLAP correction we set XAmax = 1-0 
while for the others we set XAma,x = 0.4. Further we use ^Tmin = 1 GeV for DGLAP models and 
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central jet 


central jet (scaled) 




Figure 2: The px spectra of the central jet calculated using the best fits for individual models 
versus the CMS data. For the bottom plot the cross sections have been scaled by to better see 
the differences between the models. 


^Tmin = 0.1 GeV for the others. Finally, we use fc^max = 100 GeV for exponential initial condition 
and ^Tmax = 400 GeV for the others. The last comment concerns the hard scale choice: in all fits 
we have used the average px of the jets. 

The influence of the Sudakov resummation model is illustrated in Fig.|^ Here, we have chosen 
the best fits to illustrate the effect. We see, that the jet spectra are rather weakly affected by the 
resummation, although the forward jet spectrum becomes steeper than the data. 

The obtained UGDs are plotted in one-dimensional plots in Fig. [^as a function of x and kx- 
Note, that in order to better reflect the difference between UGDs we plot T{x, kx)- We show 
results of all the models of Table hence also those with rather high value (see Table . All 
the UGDs with the DGLAP contribution are comparable, which shows that the evolution scenario 
is more important than a particular shape of the initial parametrization. On the other hand, the 
differences between UGDs are more pronounced than those in the px spectra, which means that 
the currently available data are not sufficient to discriminate among the models. The two most 
differing UGDs correspond to the BFKL+G (POW) and BFKL+GR (EXP) models which however 
have significantly higher y^NDP (above 1.5). 

We compare the new LHG-based UGDs with the one evolving according to a complicated 
evolution of |441 155] and fitted to HERA data [5| (we abbreviate it as ’KS-HERA’ on the figure). 
This evolution equation contains the kinematic constraint, full DGLAP correction (including quarks 
via coupled equations) and a nonlinear term motivated by the Balitsky-Kovchegov equation. The 
Px spectra resulting from this gluon density are presented in Fig. 


6 Azimuthal decorrelations 

In order to apply the fits in practice we have calculated another observable for central-forward dijet 
production, namely, the differential cross sections in azimuthal angle Acj) between the two jets. At 
leading order the two jets are produced exactly back-to-back and the distribution is the Dirac 
delta at Acj) = ir. However, due to QGD emissions of additional partons (either forming additional 
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Figure 3: The pt spectra of the forward jet calculated using the best fits for individual models 
versus the CMS data. For the bottom plot the cross sections have been scaled by to better see 
the differences between the models. 
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Figure 4: An effect of the Sudakov resummation model (BFKL+CDS) when applied to one of our 
fits for the model BFKL+CD with exponential initial condition. For comparison we plot also the 
spectra obtained from the unintegrated gluon density with more involved evolution and fitted to 
HERA data (KS-HERA), see the main text for more details. 
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Figure 5: Unintegrated gluon distributions evolving due to the niodels[T][^with the initial conditions 
A||C obtained from the fits to the LHC data as a function of x (top) and kx (bottom). The UGDs 
are multiplied by to better illustrate the differences between the models. The most differing 
UGDs are those without the DGLAP correction and with significantly higher x^/NDP > 1.5 
(BFKL+C and BFKL+RC). 
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jets or being soft particles with small pt) the two jets are decorrelated. On theory side these 
decorrelations are well described by QCD-based parton shower algorithms. However, within the 
HEF there is a natural decorrelation mechanism built-in. Namely, due to the internal transverse 
momentum kx of a gluon the dijet system with transverse momenta Pti, Pt2 is unbalanced by 
the amount \pti +Pti| = = kx- One can think of kx as a cumulative transverse momentum 

of many gluon emissions. In general, these emissions can be small-pT and large-pT emissions as 
well. The large-pT emissions may in general contribute a jet, thus we consider an inclusive dijet 
observables. 

Using the new fits and the Lx Jet program we have calculated the azimuthal decorrelations for 
the kinematics described in the beginning of Section The results are presented in Fig. The 
bands represent uncertainty that comes from the scale variation by a factor of two. We compare 
our calculation with the preliminary CMS data eQ 


7 Discussion 

In the present paper we went through a thorough study of various small-a; evolution equations 
analyzing an impact of various effects on jet observables. The effects we mean here, are certain sub¬ 
leading corrections to the BFKL equation, such as the kinematic constraint or DGLAP corrections. 
Our study was based on fitting these evolution scenarios to two samples of LHC data for high-pr 
spectra for dijet production. These samples consist of separate spectra for the central rapidity and 
forward rapidity jets. 

Our findings can be summarized as follows. First observation is that both forward jet and cen¬ 
tral jet spectra can be simultaneously and reasonably described by the High Energy Factorization 
approach and BFKL-like evolution. We obtain the best quality fits for BFKL with DGLAP correc¬ 
tion and kinematic constraint, with the DGLAP correction being the most important additional 
ingredient. This matches the fact that the data under consideration can be nicely described by 
the collinear factorization with a parton shower Whereas in the High Energy Factorization 

the parton shower is - to some extent - simulated by the transverse momentum dependent gluon 
distribution with the DGLAP correction. For all evolution models we get very good fits to the 
central jet spectrum, while most of the models have problems with precise reproduction of the 
shape of the forward jet spectrum. Several models properly describe the dijet data despite some 
differences in the resulting UGDs. Measurements of some other observables or more differential 
dijet data could help to discriminate among the models. 

Using our fits we have calculated azimuthal decorrelations for the same kinematic domain. 
This observable was also measured by GMS. The comparison of our calculation with the data is 
reasonably good, especially when using the Sudakov resummation model on the top of the evolution 
models. Interestingly, the same resummation procedure spoils the forward jet px spectrum. 

Our final remark is that although the High Energy Factorization with improved BFKL evolution 
equation catches the main physical aspects of the jet production at small x, one definitely needs 
higher order corrections. Such calculations exist for certain small x processes like Mueller-Navelet 
jets |451I55] or inclusive hadron production p+A collisions within GGG formalism |601 ICT] , but not 
for the high-pr dijet observables under consideration. 

^We note that the total cross section obtained from does not agree with [8]. The ratio of the two is approx. 1.8. 
If this is a normalization difference only, our predictions should be shifted up by this factor. 
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Figure 6: The results for the azimuthal decorrelations for inclusive forward-central dijet production 
using our best fits. When the Sudakov resummation model is applied to the generated events we 
get a better description of the CMS data. 
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